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MATTHIAS BECK 
Dedicated to Lou Billera on the occasion of his sixtieth birthday. 

Abstract. We present a new tool to compute the number <^A(b) of integer solutions to the linear 
system 

X > , Ax = b , 

where the coefRcients of A and b are integral. (/'A(b) is often described as a vector partition 
function. Our methods use partial fraction expansions of Euler's generating function for </)A(b). A 
special class of vector partition functions are Ehrhart (quasi-)polynomials counting integer points 
in dilated polytopes. 



1. Euler's generating function 

We are interested in computing the number of integer solutions of the linear system 

X G , Ax = b , 

where A is a nonnegative (m x (i)-integral matrix and b € Z™'. We think of A as fixed and study 
the number of solutions (I)a{)^) as a function of b. (Strictly speaking, this function should only be 
defined for those b which lie in the nonnegative linear span of the columns of A.) 

The function (pp^ih), often called a vector partition function, appears in a wealth of mathematical 
areas and beyond: Number Theory (partitions), Discrete Geometry (polyhedra), Commutative Al- 
gebra (Hilbert series), Algebraic Geometry (toric varieties), Representation Theory (tensor product 
multiplicities), Optimization (integer programming), as well as applications to Chemistry, Biology, 
Physics, Computer Science, and Economics. 

Denote the columns of A by ci, . . . , c^. The following lemma goes back to at least Euler |15j : 
Lemma 1 (Euler). </>A(b) equals the coefficient of := z^^ ■ ■ ■ of the function 

■^^^^ " (l-Z^l)---(l-Z^d) 
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expanded as a power series centered at z = 0. 



Proof. Expand each factor of the right-hand side into a geometric series. □ 

Equivalently, the coefficient of in /(z) equals the constant term in "^e^, denoted by const 
So Euler's Lemma can be conveniently stated as 

'^^^^^ = ™* (i-zco...a-zc.)zb ■ 

In a series of articles jSl El IHl Ej , we used complex integration of to compute 4>a (b) for special 

cases of A. Similar techniques were applied in jl9j . Here we expand into partial fractions to 
compute its constant term, and hence (j)AO^)- This work constitutes in a sense a refinement of the 
complex-integration methods, with the advantage that it is more flexible and — more importantly — 
applicable for any integral linear system. 



2. Vector partition functions 



The nature of the counting function (/)A(b) is given by the following theorem. A quasi-polynomial 
is a finite sum Q{h) = X^nCn(b) b" with coefficients Cn that are functions of b which are periodic 
in every component of b. The degree of Q is the degree of the largest power b" appearing in Q. A 
matrix is unimodular if every square submatrix has determinant ±1. 

Theorem 2 (Sturmfels [22] )• The function (/'A(b) is a piecewise- defined quasi-polynomial in b of 
degree d — rank(A). The regions o/M™ in which 0a (b) is a single quasi-polynomial are polyhedral, 
that is, they are defined by linear constraints. If A is unimodular then (j)js^ is a piecewise- defined 
polynomial. 



The unimodular case of this theorem is due to Dahmen and Micchelli III] Corollary 3.1]. 

The computation of both i;^A(b) and the chamber complex consisting of the regions of quasi- 
polynomiality give rise to challenging problems. The most powerful technique for computing 0a (b) 
which we are aware of is due to Brion and Vergne [S]. (The methods described here are much more 
elementary.) The chamber complex is still much of a mystery. A promising approach can be found 
in HSI- 

We finish this section with a reciprocity theorem. Let (/>^(b) count the integer solutions of 

X > , Ax = b . 

Both and are quasi-polynomials, and can hence be algebraically defined for arguments which 
are not integer vectors in the positive span of A. The following identity shows the close relationship 
between the two functions. 

Theorem 3(4). The quasi-polynomials 0a (^nd satisfy 

0A(-b) = (-l)'^--°'^^0l(b) . 
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This identity gives rise to a symmetry property of (/>A(b). Let denote the sum of the entries in 
the k^^ row of A, and let r = (ri, . . . , r^). Then the integer solutions of 

X > , Ax = b 

are in bijection (via ^ — 1) with the integer solutions of 

x>0, Ax = b — r, 

and hence (/>^(b) = </)A(b — r). This yields: 

Corollary 4. The quasi-polynomial (pA. satisfies 

<^A(b) = (-l)'^-'^^°'^^</.A(-b-r) . 



3. The partial-fractions method 



This section describes the idea behind our computations. Recall that our goal is to derive 

4>a(^) = const r ^ — ^ . 

' (1-Z^l)---(1-Z^d)zb 

We start by expanding 

1 



(1 - Z^l) • • • (1 - Z''d)z^ 

into partial fractions in one of the components of z, say zi: 



1 J^Ak{z,h) , Xb,{z) 



(l_zci)...(l_zcd)zb l_zCfe /L. 

Here and Bj are polynomials in zi, rational functions in Z2, Zm, and exponential in bi. (We 
are tacitly assuming that there are no multiple poles (besides z = 0), which generally holds unless 
m = 1, a case which can be handled easily.) The two sums on the right-hand side correspond to 
the analytic and the meromorphic part with respect to zi = 0. The latter does not contribute to 
the zi-constant term, whence 

1 , fJ^Akiz,h)\\ 



6a (b) = const^^,...,;,^ — — const^i V 



vfc=l 



(1) = const 



The effect of one partial fraction expansion is to eliminate one of the variables of the generating 
function, at the cost of replacing one rational function by a sum of such. The following idea is now 
evident. 

Algorithm. Apply ^ repeatedly to eliminate Z\, then Z2, etc., up to z^-i- 



The constant term of the remaining rational functions in the one variable can be computed with 
the methods introduced in 5 . This implies, in particular, that any vector partition function is a 
quasi-polynomial whose nontrivial ingredients are Fourier- Dedekind sums 5 . The constant-term 
computation for the last variable also explains the quasi-polynomial character of (j)A, since in this 
last step roots of unity will appear, with the components of b as exponents. 
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We note that this algorithm computes 0a (b) as a function ofh, that is, it allows symbolic compu- 
tation. Secondly, it is not very hard to deduce Sturmfels's Theorem |2l from this algorithm. What 
might be more important, however, is the fact that the constraints which define the regions of 
quasi-polynomiality of (pA are obtained "on the go" as one computes 4>a- When expanding into 
partial fractions, one has to check where the poles of a rational function are. The components of b 
will appear (linearly) in the exponents of these rational functions, and hence one will automatically 
have to split the computation into cases which give rise to different regions in the chamber complex 
of quasi-polynomiality. 

While we hope that our algorithm gives a computational tool, in particular, for obtaining the 
chamber complex defining the regions of quasi-polynomiality of a vector partition function, it is not 
clear to us why the chambers obtained from the partial fractions analysis coincide with Sturmfels's 
predicted chambers. 

Our algorithm is best illustrated by going through an actual example. 



4. An illustrating example 



Let A=^| ^ Q l)' '^rite b = (a, b), so 0a(«> b) counts the integer solutions of 



xi + 2x2 + X3 = a 



By Euler's Lemma ^ 

1 



(t>A(cL, b) = const r-. r, N , N / N !. • 

' ^ {1 - zw){l - z^w){l - z){l - w)z''w'' 

We first expand into partial fractions with respect to w: 



1 (1-Z)2 (l-2)(l-22) 

T 

k=l 



{\ — zvS)(\ — z'^w)(\ — vS)w^ 1 — ZW 1 — z'^w 1 — W j-^^ 

Taking constant terms gives 



4>A{ci,b) = constj; ( — — const^ 



(1 — z)z'^ \(1 ~ zw){l — z'^w){l — w)w^ 
const — -TT -I- — — TT- + 



^26— a+3 

(2) = const ( -— -17 -I- 7" --Kj- ^ -I- 



(1-Z)3 (1_2)2(1_^2) (l-z)2(l_22)^a 

(At this point, we could interpret each of the three constant terms as counting integer solutions 
to new linear systems. For example, the last term gives the number of integer points (x, y) > 
satisfying 2x + y < a. To keep a general flavor, we continue with our general algorithm.) 
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We compute the constant term of each of the three terms separately. For starters, 

^b—a+l 

const -. -5- = 

if & — a + 1 > 0, equivalently (since a and b are integers) b > a. If 6 < a, we use 



k + 2\ 1^ 



which gives 



(l-z)3 2 



1 fa-b + l\ ia-bf a-b 
const — T— r = „ = h 

(1 - z)3z«-^-i \ 2 



Incidentally, this identity is true not only for b < a, but also for b = a and b = a + 1, because 
the right-hand side vanishes then. This suggests that the regions of quasi-polynomiality overlap, 
as was proved by Szenes and Vergne Hence for the first constant term, we obtain 

const 



1 if6>a. 



(l_^)3^a-fe-l ]i^^a^ if6<a + l. 



For the second term in (jSJ, 



^2b-a+3 



if 26 — a + 2 > 0. If a > 26 + 3 we expand into partial fractions again: 

const — — = const ' + — ^ -rr^ + ^ ^ ^ +^ ' ' 



(l-z)2(l-z2) \^(1_^)3 ' (1_^)2 ' Y-Z 1 + Z J 

_ (a -26)2 26 -a 1 + (-1)"+^ 
~ 4 ^ 2 ^ 8 ■ 

Similar to the first constant-term computation, this identity is also valid for a = 26 -|- 2, 26 -|- 1, 26, 
whence 

z2fe-«+3 ^ Jq if a < 26 2, 

The computation for the last term in @ is almost identical. (Note that this term always con- 
tributes.) 

1 /l/2 f + i T + f + l (-1)V8\ 
const -TTTT rr^ — = const ' + ^ + -^—^ + ^ ' ' ' 



(1 - Z)2(l - z2)2« 1(1-2)3 (l-z)2 l-z l + Z 
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Summing up all terms in gives finally 
4>A{a,b) -- 



4 -t- U -I- g 



2 



o ~r /I 



+ 



2b-a 



4 

a6 



g 
4 



2 



^ + ^ + 1 



a+b , T+C-l)" 
2 8 



+ 



7+(-l)° 



2 ' 8 

if a < 6, 

iff-l<6<a + l, 



if 6 < f. 



if a < 6, 

iff-l<6<a + l, 
if5<f 



It is a fun exercise to show that </>a(o, b) = (/'a(— « — 4, —6 — 3), as promised by Corollary^] 

We pause for a moment to introduce the geometry behind our linear system. By thinking of X3 
and 2:4 as slack variables, we can see that 



Xi + 2X2 + X3 
Xl+ X2 + X4 



a 
b . 



is equivalent to 



Xi,X2 > , 



Xi + 2X2 < 
Xi+ X2 < b . 



Depending on the relationship between a and b, the geometric figure described by this linear system 
is a quadrilateral or triangle. 



(0,a/2) 



(0,0) 



(26 -a, a -6) 




The inequalities which define the different cases for (j)p^{a,b) determine the chamber complex in 
the parameter space, and in this example also the combinatorial type of the polytope: Except for 
boundary cases, in the first (a < b) and last {b < |) case, it is a triangle, and in the second case 
(a > 6 > |) a quadrilateral. In fact, the geometry of the triangle in the first case only depends 
on a, which is reflected in (/>a(oj^)j and similarly the last case shows only dependency on b. In 
this last case, the triangle actually has integer vertices (the linear system is unimodular), whence 
we get a polynomial. In the first two cases the vertices of the polytope are half-integral, which is 
reflected in the period-2 quasi-polynomials. The following picture illustrates the three overlapping 
chambers. 




5. Ehrhart quasi-polynomials 

A convex polytope V in R'^ is the convex hull of finitely many points in R"'. Alternatively (and this 
correspondence is nontrivial P^)' define V as the bounded intersection of affine halfspaces. 

A polytope is rational if all of its vertices have rational coordinates. (A vertex of P is a point v G P 
for which there is a hyperplane H such that {f } = V H.) We denote by V° the relative interior 
of V. For a positive integer t, let L-p^t) denote the number of integer points ("lattice points") in 
the dilated polytope tV = {tx : x £ V}. The fundamental result about the structure of L-p is as 
follows. 

Theorem 5 (Ehrhart Jl]). If V is a convex rational polytope, then the functions L-p{t) and L-p°{t) 
are quasi-polynomials in t whose degree is the dimension ofV. If V has integer vertices, then L-p 
and Lpo are polynomials. 

Ehrhart conjectured and partially proved the following reciprocity law, which was proved by Mac- 
donald HO]. 

Theorem 6 (Ehrhart-Macdonald) . The quasi-polynomials Lp and Lp° satisfy 

Lp{-t) = {-lf^'PLp.{t) . 

The computation of Ehrhart quasi-polynomials is only slightly easier than that of (p^. ■ Recent work 

includes [niHiisiiniiiiiiiiiiniiiziiiaiii]- 

Suppose the convex polytope P C R'^ is given by an intersection of halfspaces, that is, 

= |x G R"' : Ax < b| , 

for some (m x d)-matrix A and m-dimensional vector b. We may convert these inequalities into 
equalities by introducing slack variables. If V has rational vertices, we can choose A and b in such 
a way that all their entries are integers, without loss of generality nonnegative ones. In summary, 
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we may assume that a convex rational polytope V is given by 

7' = {xGM^o: Ax = b} , 

where A € Mmxd(^>o) ^-nd b € Z>o- (If we are interested in counting the integer points in V, we 
may assume that V is in the nonnegative orthant, i.e., the points in V have nonnegative coordinates, 
as translation by an integer vector does not change the lattice-point count.) 

The connection to vector partition functions is now evident. Since tV = {x € M^q : Ax = tb}, 
we obtain L-p{t) = (f)Aith) as a special evaluation of (pA_. Note that th,t = 1,2,... lie in the 
same chamber of quasi-polynomiality of 0a- Ehrhart's Theorem El is therefore a special case of 
Sturmfels's Theorem [21 and Theorem (HI is a special case of Theorem |S1 

As an example, the quadrilateral Q described by 

x + 2y < 5 
x,y> , Z- A 

(a special case of the polygons appearing in Section^ with vertices (0, 0), (4, 0), (3, 1), (0, 5/2) has 
the Ehrhart quasi-polynomial 

Ls(t) = c/,A(5t, 4t) = ^ t2 ^ ? t + '^+1"^^' . 



Ehrhart quasi-polynomials are easier to compute, since one does not need to derive/know the 
chamber complex of quasi-polynomiality of ^p^. Our algorithm, naturally, works just as well for 



6. Concluding remarks 



Many open problems and questions remain untouched. Most importantly, we hope that our ideas 
will give rise to practical implementations of computing vector partition functions (including the 
chamber complex defining the regions of quasi-polynomiality) in their various disguises in form of 
a computer software. Again, the partial-fraction method is more general and more flexible than 
algorithms using complex residues, for example in 16 . 

Not unrelated is the question of computational complexity. It is known that the rational generating 
function of the Ehrhart quasi-polynomial of a d-polytope can be computed in polynomial time if d 
is fixed |21. We have not analyzed the complexity of our algorithm. This seems interesting in the 
light that our algorithm depends on the number of contraints, geometrically corresponding to the 
facets (codimension-1 faces) of the polytope. Barvinok's algorithm, in contrast, depends on the 
number of vertices of the polytope. 

Another venue which should be explored is the following: After each step in our algorithm, one 
could reinterpret the constant term of each summand as counting integer solutions of a linear 
system. It might be interesting to simplify each of these linear systems (as far as this is possible 
without changing the number of integer solutions), and then proceed with the algorithm. 
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